A new phase of scalar field in a non-minimal model for inflation 
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W consider the dynamics of a scalar field non-minimally coupled to gravity in the context of 
cosmology. It is demonstrated that there exists a new phase for the scalar field, in addition to the 
infiationary and dust-like (reheating period) phases. Analytic expressions for the scalar field and 
the Hubble parameter, which describe the new phase are given. The Hubble parameter indicates 
an accelerating expanding Universe. We explicitly show that the scalar field oscillates with time- 
dependent frequency. Moreover, an interaction between the scalar field in the new phase and other 
fields is discussed. It turns out that the parametric resonance is absent, which is another crucial 
difference between the dynamic of the scalar field in the new phase and dust-like phase. 



PACS numbers: 98.80.Cq 



I. INTRODUCTION 



To solve the flatness problem and the horizon prob- 
lem in cosmology, Alan Guth introduced the inflation 
paradigm The simple model, which is a single scalar 
field (inflaton) with minimally coupling to gravity, is de- 
scribed by an action 



5* 



, (1) 



where = SttG. Although model ([T]) with quadratic 
potential, V{tp) = l/2m^(^^, is consistent with the cos- 
mic microwave background [4j, many models have been 
proposed to produce the inflation era . Regarding the 
motivations, many attempts are devoted to derive the 
inflation era by more fundamental principles or find con- 
nections between inflaton and other fields that are used 
in particle physics(e.g. Higgs field) [1, H, Q. 
An intuitive picture of the dynamics of inflaton is as fol- 
lows: the inflaton field rolls slowly down its potential (in- 
flationary period), eventually the field(s)oscillates around 
the minimum of its potential and decays into light par- 
ticles (reheating period). It has been assumed that the 
reheating period takes place just after the inflationary 
period. Also, since the inflaton has lost its energy, it be- 
haves like dust matter in the reheating period. 
For the quadratic potential, this intuitive picture is sup- 
ported by analytic methods ( see Q for the quartic 
potential). Usually model builders use the stated sce- 
nario for their models ( see Q for a different scenario). 
Then, they set limits on the parameters of the models. 
So, one can divide the dynamic of inflaton in two phases; 
inflationary phase and matter dominate phase(near min- 
imum of the potential) . 

In this work, we show that in the following effective ac- 
tion 



S 



(2) 

where G'^'^ is the Einstein's tensor, the scalar field has 
a new phase. Germani et.al Q, proposed the above ac- 
tion with quartic potential, A(p^, in the context of the 



new Higgs inflation. The cosmological perturbation of 
the model and reheating period of the model have been 
studied in [l^l and [ll| respectively. 
This paper is organized as follows: in §11 we briefly re- 
view the model and obtain equations, then we qualita- 
tively discuss why the scalar field of model has a 
new phase for any typical potential, I^(</3). In §111 ana- 
lytic solutions for the Hubble parameter and the scalar 
field are provided for the quadratic potential. The solu- 
tions describe the new phase of the scalar field. In §IV 
we investigate an implication of the new phase by con- 
sidering an interaction between a relativistic field and 
the scalar field. To introduce reader to the method that 
is used in this paper, we re-derive solution of the scalar 
field of the model ([T]) in reheating period in Appendix A, 
which is identical with the solution in the textbook [5|. 
Appendix B is devoted to some properties of the Fresnel 
integrals. 



II. EVIDENCES FOR EXISTENCE OF THE 
NEW PHASE 

To obtain field equation and the Friedman equation of 
the model in FRW background metric, we use the ADM 
formalism with the following metric ansatz [l^ 

ds^ = -N{t)dt^ + a{tf5ijdx'dx^ . (3) 

By inserting this expression into ^ we have [lo| 

ij2 1(^2 SH^a'^if'^ 
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(4) 

by varying the actionQ with respect to the laps N and 
ip and setting A'^ to 1, we obtain 

,.2 



dip 
(5) 



The behavior of the scalar field can be divided into two 
regimes 



2 



If aH ^ 1, from (O we have 



ij2 
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(p + iHip = — 



[^2 + 2V{^)\ 



dip 



(6a) 
(6b) 



Equations ()6ap and (j6bp are the same as equations which 
arc derived from ([T]) and studied in textbooks [IjQ- Since 
we want to compare other regime with this regime, we 
quote the main results. 

The second term on the left hand side of (|6b[) . which is 
always positive, acts as a dissipativc force. Therefore, 
the scalar filed rolls toward the minimum of the poten- 
tial, and it behaves as a dust matter. For the simplest 
potential, y(iy9) = l/2m^iy9^, it has been shown that 



Hit) = — 
^ ' 'it 



cpit) 



cos [mt + 5(0)] (7) 



where g{0) is an arbitrarily constant. It is important to 
note that expressions in (O arc only valid for H{t) < m, 
so higher order terms, have been neglected [Bj. 



• If aH ^ 1, from ([S]) we have 



H' 
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[(p^9a^H'^ 



ip — SHipw 
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2V{ip)] , 



3H^a^ dip ' 
where if is the effective equation of state 
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(8b) 



(9) 



The inflationary phase of this regime, u; — > — 1, was stud- 
ied in 0, [l^. Note that during the inflationary phase, 
the second term on the left hand side in (j8b[) is posi- 
tive, and acts as a dissipative force. So, again, the scalar 
filed rolls toward the minimum of the potential , but at 
H = 2/3t, it is vanished and, after some time, its sign 
may be changed and act as a driving force. Another point 
is that the right hand side of (|8b[) . depends on H , which 
is not constant in this period. This behavior, shows that 
we need a different analysis for this regime. 



III. EXPLICIT SOLUTIONS FOR THE 
QUADRATIC POTENTIAL 

In this section we consider the simple potential V{(p) = 
\/2rn?ip'^, then we seek solutions for the equations. We 
use the so-called averaging method [31 , which is used in 
study of nonlinear differential equations. In Appendix A, 
we use this method to re-derive (O. 
Let us define a new time variable t as 



t = 



3aH 



dr. 



(10) 



It is worth to mention that this step is only a trick to 
solve the equations, when we obtain the solutions we'll 
back to the original time variable, t. 
Combining Eqs. ([5a|) . (j8bp and (|10l) . with the quadratic 
potential, gives 



H = — m (^ 



H m 



(11) 



-3(/), 



where prime denotes the differential with respect to r. 
In order to use the method of averaging, we define the 
following relations [l^ 



/ 



cos (r -t- /(r)) ^ 

Km 

-V6H 



■ sin (t + /(t)) 



(12a) 
(12b) 



where /(t) is an arbitrary function. Differentiation of the 
righthand side of (|12ap must be equals to the righthand 
side of (jl2bp . Hence we have 

H' cos (r + fir)) = Hf sin (r + /(r)) , (13) 

Substitution of (|12ap and (jl2bp into ([TT]) results in 



2H' sin (r -t- /(r)) + Hf cos (r + /(t)) 

+ sin (r + /(r)) - 2H cos (r + /(r)) = 0. 



(14) 



H' and /' can be found from Equations and ([TH]) as 



H' 



cos (T + /(r)) 
1 + sin^ (r + /(r)) 

g sin (t + /(t)) 
' l + sin^T + fir)) 



2 cos (r + /(T)) 



2 cos (r + /(r)) 
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sin(T-H/(r)) 



■sin (t + /(t)) 



(15) 



So far we have not used any approximation for and 
/(r). In the method of averaging, equations in are 
replaced by averaged expressions. For this goal, note that 
if we keep H and / fixed, the righthand side of (fTS)) are 
TT-periodic in r, so we can average over r. 
By applying the stated procedure, one can show that 



\/2- 1, 



cos^ (r + /(t)) 
l + sin2 (t + /(t)) 
cos (t + /(t)) sin (t + f{T)) 
l + sin2 (r + /(T)) 



sin^ (r + /(r)) 
1 + sin2 (r + /(r)) 



2 ' 



0, 



(16) 



where < ••• >= <j> {■ ■ ■ )dT / it denotes an average over r 
(but H and / arc fixed). So, from Eqs. and (|16p . 
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we have 



IV. AN IMPLICATION 



/' = 2(x/2-l) 
H' - 



9H^a V2 
m 2 



(17) 



The averaged equations are very easy to solve: 



In this section, to show an imphcation of properties of 
the new phase, we wih consider an interaction between a 
matter field, x, and the scalar field. 

The decay of the scalar field into a relativistic field can 
be described by ^ 



I8a. V2, 
(1 - -TtF 

771 2 



/(r)- 2(^/2- l)r. (18) 



Hir) 



We can use formula ([10)) and p8)) to obtain explicit re- 
lation between t and t as 



t 



a 2 



r2 . 



By using Eqs. (fTS)) and (|19p we have 

H{t) = 

3(2 - V2)t 

so, the effective equation of state becomes 

2 H 

w = -l ^ w -0.41. 

3^2 

Moreover, from (|12ap and ([TO]) , we also obtain 



, , V6H(t) 
ip{t) = cos 



£(2-V2)(y2-^)t2 



(19) 



(20) 



(21) 



(22) 



where H{t) is given by ([20)1 . 

The new phase of the scalar field is described by ([^0)1 and 
(|22|) . According to ([21]) , for this new phase we have w < 
— 1/3, so, the Universe is driven to accelerated expansion 
phase by the scalar field. Expression ([^^ . shows that the 
scalar field oscillates with time-dependent " frequency" , 
as indicated in Fig.l. Both of the properties are more 
different than the dust-like phase, which is described by 



^222 

2^ V X ■ 



(23) 



It has been shown that if the interaction takes place dur- 
ing the dust-like phase, the parametric resonance insta- 
bility is occurred d, [l^ ■ 

Notice that this result is obtained if the expansion of the 
universe is neglected. The stated assumption, seems to 
be a reasonable condition, at least for the first approx- 
imation, if the rate of interaction is too fast compared 
to the Hubble expansion time, see discussions about this 
note in Q. As mentioned in [Ej, what one actually finds 
from present of the parametric resonance, is that the per- 
turbative analysis is rather misleading, during dust-like 
phase of the infiaton. 

Here, our aim is to quest for the parametric resonance, 
during the new phase of the scalar field. 
If the relativistic field, Xj is decomposed into Fourier 
modes as 



{t)e 



-ik.X 



(24) 

then using ([22]) and ([23)) . the following equation for the 
Fourier modes is obtained 



Xk + (fc' + g'^A^ cos^iMh^)) xk - 0, (25) 



where 



.9 6 i?^ „ 



,2^2' M^^^(2-V2)(x/2--). (26) 



Now, following Q , we will neglect the expansion of space. 
So, If we define the following variables 




FIG. 1. The scalar field versus t for nm — 10, m — 10a. 



1,2 , 2 a2 2/12 



2AP 



4AP 



it thus follows that 
d^Xk 



+ (a;2 + £cos^2)xfe = 0. 



(28) 



To obtain an approximate solution for Xk , we expand Xk 
as 



Xk = Xk + £xl + 



(29) 



The expression ([29)) is valid, if Xk has not terms that grow 
without bound as t — > oo. So, it is necessary to cross- 
check this condition after we obtain an explicit expression 

for xl- 



4 



Substituting ((29)) into (|28p . and keeping all terms up to 
second-order in e, yields 



. , ,2 
— + ^ Xk 







-Xlcose 



(30) 



The first equation in (pO| gives 

6i sin uj^ + 62 cos w^, 



(31) 



where 61 and &2 are the integration constants. 
Using (pij) and the following variables 



(32) 



the second equation in ((5D|) can be solved as 

xl =Qi cos(w^) + Q2 sin(a;^) - -^^j^ 

cos(w^ + a;2)[6^C(u) - b2S{u)] + 
sm{oj£, + oj'^)[b,C{u) + biS{u)] + 
cos{uj^ - uj^)[b,C{v) + b2S{v)] + 
sm{u^ - oj^)[b,C{v) - biS{v)] + 



(33) 



2C{- 



-)[b2 sin(a;^) — bi cos(aj^)] 



Where Qi, Q2, are the integration constants, and C{x), 
S{x) are the Fresnel integrals ( see Appendix B). 
From properties of the Fresnel integrals [l5|, we conclude 
that expression p3p has not any terms that grow without 
bound as i — > 00. Therefore, we expect that the pertur- 
bation approximation for x, which is given by (|29p . is 
assured for all times. 

Figure 2 compares the approximate solution for x, to the 
numerical solution. The two curves are almost indistin- 
guishable. 

If we used ([25)) instead of (j6bp , we would have terms that 
grow without bound as i — ^ 00, which are sources of the 
parametric resonance during dust-like phase of the scalar 
field. 
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Appendix A 

In this appendix, to show the power of the method that 
is used in this paper, we re-derive ([7]), by the averaging 
method. 




(b) 



FIG. 2. The relativistic field, x, versus ^. Solid line represents 
the numerical solution. Dashed line represents the approxi- 
mate solution. Where the parameters are chosen such that 

X = 16.5, — = 12. The two curves are so similar that it is 

alm ost impossible to distinguish between them, as shown in 

m 



By introducing /3 = mt, and V{lp) = l/2'm?tf'^ , equations 
([6a| and ()6bp can be written to give 



H 



(Al) 



(p + 3 Lp = — (yS . 

m 

Now, consider the following relations 

^ = ^cos(/3 + 5(/3)), 
Kin 

¥'' = ^^sin(/3 + .9(/3)). 



(A2a) 
(A2b) 



where g{P) is an arbitrary function. Differentiation of the 
righthand side of ()A2ap must be equals to the righthand 
side of (|A2bp . Hence we have 



H' cos {P -f g(fi)) = Hg' sin (/3 + 5(/?)) , 



(A3) 
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Substitution of (|A2a|) and (IA2bp into (KA) yield 

H' sin (/3 + + Hg' cos (/3 + g(/3)) 



+ — sin(/3 + 5(/3)) = 0. 



(A4) 



By algebraic manipulations, H' and can be found from 
Equations (K3\i and ((X4| as 



.9 = - ^ sin (/3 + 5(/3)) cos (/3 + 5(/3)) 



(A5) 



sin^ (/? + .9(/3)). 



These expression for H, and are exact. 
The advantage of the averaging method is that an ap- 
proximation to the solution of (|A5|) can be obtained by 
replacing (|A5p with its averaged equations as follows: if 
we keeping H and g fixed, the righthand side of (|A5p are 
TT-periodic in /3. By noting that 

(sin^ (/3 + g(/3))) = 1 (cos (13 + sin (/? + giP))} 



where H(t) is given by (|A7p . 

According to the averaging method, the above expres- 
sions for Hubble parameter and the scalar field are valid 
for H{t) < m. The results are the same as those obtained 
in 0, which are also valid for H{t) < m. 



Appendix B 

The Fresnel integrals are defined by [lE 

C{x) = / cos(-7rx^)dx 
Jo 2 

1 

S{x) = / sin(-7ra; )dx 
Jo 2 

A series expansion for a; < 1 gives 
— (-l)"(7r/2)2" 



(Bl) 



(A6) 

where <> is used to indicate that we average over /?, the 
averaged equations can be solved as 



H{t) = Yt ' 9{t)=9{0), 



(A7) 



n (B2) 

^ (2n+l)!(4n + 3) 

n— 

Therefore, lim C(x) = lim S{x) = 0. 
Asymptotic expansion of the integral are given by 

1 1 . ,1 



where g(0) is an arbitrary constant. Using (jA2ap and 
(|A7p . we have 



Cix) = ^ + — sin(^7ra;2), 

2 TTX 2 

S'(x) = i - — cos(i7rx2). 
2 7r.T 2 



Hence, lim C{x) ^ lim S{x) = 1/2. 



(B3) 
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